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Local  atomic  structure  of  Ca-Mg-Zn  metallic  glasses 
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The  amorphous  structure  of  four  Ca60MgxZn40_x  (X=10,  15,  20,  and  25  at.  %)  ternary  metallic  glasses 
(MGs)  has  been  investigated  by  neutron  and  x-ray  diffraction,  using  Reverse  Monte  Carlo  modeling  to  simu¬ 
late  the  results.  A  critical  analysis  of  the  resultant  models,  corroborated  by  ab  initio  molecular-dynamics 
simulations,  indicate  that  the  glass  structure  for  this  system  can  be  described  as  a  mixture  of  Mg-  and  Zn- 
centered  clusters,  with  Ca  dominating  in  the  first  coordination  shell  of  these  clusters.  A  coordination  number 
(CN)  of  10  [with  about  7  Ca  and  3  (Mg+Zn)  atoms]  is  most  common  for  the  Zn-centered  clusters.  CN=  1 1  and 
12  [with  about  7-8  Ca  and  4  (Mg+Zn)  atoms]  are  most  common  for  Mg-centered  clusters.  Fivefold  bond 
configurations  (pentagonal  pyramids)  dominate  (—60%)  in  the  first  coordination  shell  of  the  clusters,  suggest¬ 
ing  dense  atomic  packing.  Bond-angle  distributions  suggest  near-equilateral  triangles  and  pentagonal  bipyra¬ 
mids  to  be  the  most  common  nearest  atom  configurations.  This  is  the  systematic  characterization  of  the 
structure  of  Ca-Mg-Zn  MGs,  a  category  of  bulk  MGs  with  interesting  properties  and  intriguing  applications.  It 
is  also  the  experimental  verification  of  the  principle  of  efficient  packing  of  solute-centered  clusters  in  ternary 
MGs. 

DOI:  10.11 03/Phy sRevB .82.1 04206  PACS  number(s):  61.43.Dq,  61.05.fm,  61.43.Bn,  83.10.Rs 


I.  INTRODUCTION 

The  recently  discovered  ternary  Ca-Mg-Zn,  Ca-Mg-Cu, 
and  Ca-Mg-Al  bulk  metallic  glasses  (BMGs)  have  unique 
properties  that  distinguish  them  from  other  BMGs.1-5  They 
are  based  on  two  simple  metal  elements,  Ca  and  Mg,  and 
have  the  lowest  density  (in  the  range  of  1. 6-2.4  g/cm3) 
among  all  currently  known  BMGs.6  These  metallic  glasses 
also  have  extremely  low  Young’s  (—20-30  GPa)  and  shear 
(—9-14  GPa)  moduli,  comparable  with  the  moduli  of  hu¬ 
man  bones.7  Some  of  these  ternary  alloys  have  extremely 
good  glass-forming  ability  (GFA),  similar  to  the  GFA  of  the 
best  Zr-based  BMGs.4’6,8 

Figure  1  shows  the  liquidus  projection  in  the  Ca-rich  cor¬ 
ner  of  the  Ca-Mg-Zn  system.9  The  compositions  of  the  ter¬ 
nary  Ca-Mg-Zn  amorphous  alloys  reported  in  Refs.  2  and  4 
are  marked  as  solid  circles  and  their  critical  plate  thicknesses 
are  shown  by  corresponding  numbers  (in  mm).  The  presence 
of  a  ternary  eutectic  reaction  provides  a  strong  liquidus  tem¬ 
perature  gradient  in  the  selected  composition  area,  which 
leads  to  high  sensitivity  of  GFA  on  alloy  composition.  For 
example,  the  maximum  plate  thickness  at  which  an  alloy  is 
fully  amorphous  after  casting  in  a  water-cooled  copper  mold, 
can  be  increased  from  0.5  mm  to  6-10  mm  within  a  compo¬ 
sition  range  of  only  5%. 4  Thermodynamic  analysis  of  the 
onset  driving  force  (ODF)  for  crystallization  of  different 
phases  in  the  Ca-Mg-Zn  system  has  shown10  that  the  GFA  of 
the  ternary  alloys  improves  with  a  decrease  in  the  ODF  of 
competing  intermetallic  phases  [for  example,  CaZn  (at  X 
<17.5  at.  %)  and  CaMg2  (at  Y>  17.5  at.  %)  for 
Ca60MgxZn40_x  alloys,  see  Fig.  2]. 

Recent  structural  analysis  of  metal-metal-type  BMGs  has 
focused  on  those  based  on  transition  metals  such  as  Zr-Cu 
(Refs.  11  and  12)  and  Zr-Cu- Al. 13  Many  papers  discuss  the 
icosahedral  local  structures  in  these  BMGs  and  the  impor¬ 


tance  of  these  motifs  on  relaxation  dynamics,  GFA  and  de¬ 
formation  properties.  However,  icosahedral  order  is  not  uni¬ 
versally  the  dominant  short-range  order  (SRO)  in  metallic 
glasses.  For  example,  transition-metal-metalloid  metallic 
glasses  [such  as  those  based  on  Fe-C  or  Ni-B  (Refs.  14  and 
15)],  and  many  Al-transition-metal  glasses  (with  high  Al 
content16),  are  not  characterized  by  icosahedral  SRO. 

Ca-Mg-Zn  represents  a  different  group  of  metal-metal- 
type  BMGs  based  on  alkaline-earth  metals,  for  which  the 
glass  structure  has  not  been  systematically  studied  so  far, 
although  a  few  examples  for  Mg-Cu-based  MGs  have  been 
reported.17,18  Based  on  relative  atomic  sizes19  and  the  effi- 

Mg 


FIG.  1.  Liquidus  projection  of  the  Ca-Mg-Zn  ternary  system. 
The  projection  shows  thick  solid  lines  bounding  the  Ca,  CaMg2, 
CaZn,  and  Ca3Zn  phase  fields.  Ternary  eutectic  and  peritectic  com¬ 
positions  are  marked  by  an  open  triangle  and  an  open  circle,  respec¬ 
tively.  The  compositions  of  amorphous  alloys  and  their  maximum 
plate  thicknesses  (in  mm)  are  shown  as  solid  circles  and  pertinent 
numbers.  The  trapezoid  with  a  thick  dashed  boundary  represents  the 
composition  range  for  the  formation  of  ternary  Ca-Mg-Zn  BMGs 
(Ref.  4). 
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Mole  fraction.  Mg 

FIG.  2.  Calculated  onset  driving  forces  of  various  crystalline 
phases  for  Ca60MgxZn40_x  alloys  versus  Mg  content  at  390  K.  Two 
lines  corresponding  to  maximum  onset  driving  forces  belong  to 
CaZn  and  CaMg2  phases.  The  circles  represent  the  amorphous  alloy 
compositions  with  indicated  critical  thicknesses  (Ref.  10). 

cient  solute-centered  cluster  packing  (ECP)  model  of  amor¬ 
phous  structure,20  the  Ca-Mg-Zn  ternary  glasses  are  repre¬ 
sented  as  (10,9)  efficient  cluster  packing  structures,  where 
Mg  solute  atoms  are  surrounded  by  approximately  ten  near¬ 
est  neighbors  and  Zn  solute  atoms  have  approximately  nine 
nearest  neighbors,  assuming  that  all  the  nearest  neighbors  are 
solvent  atoms  (i.e.,  Ca).  Here,  the  term  solvent  is  used  to 
indicate  the  atomic  species  in  the  glass  with  the  largest  atom 
fraction  while  the  term  solute  denotes  the  other  atomic  spe¬ 
cies  in  the  glass.  Glass  compositions  predicted  from  the  ECP 
model  range  from  Ca70Mg15Zn15  to  Ca53Mg12Zn35,21  which 
fit  within  the  trapezoid  shown  in  Fig.  1.  Structural  insights 
from  the  ECP  model  have  been  validated  in  binary  metallic 
glasses22  but  have  not  yet  been  compared  with  experimental 
data  for  ternary  metallic  glass  structures. 

It  is  the  objective  of  the  present  work  to  experimentally 
characterize  the  local  atomic  structures  in  Ca-Mg-Zn  ternary 
metallic  glasses.  In  the  following,  we  use  neutron  and  x-ray 
diffraction  data  as  the  basis  on  which  to  model  the  amor¬ 
phous  structure  of  a  series  of  ternary  Ca60MgxZn40_x  ( X 
=  10,15,20,25  at.  %)  metallic  glasses  using  Reverse  Monte 
Carlo  (RMC)  methods.  The  effect  of  alloy  composition  on 
the  structure  and  GFA  is  examined  and  the  average  amor¬ 
phous  structure  identified  in  this  work  is  compared  with  the 
amorphous  structure  predicted  by  the  ECP  model20  and  dis¬ 
cussed  in  comparison  with  those  of  other  BMG  families. 

II.  EXPERIMENTAL  PROCEDURES 
A.  Sample  preparation 

Four  amorphous  alloys,  Ca60Mg10Zn30,  Ca60Mg15Zn25, 
Ca60Mg20Zn20,  and  Ca60Mg25Zn15,  were  used  in  this  study. 
Among  these,  the  Ca60Mg10Zn30  and  Ca60Mg25Zn15  alloys 
show  marginal  GFA,  with  maximum  fully  amorphous  plate 
thicknesses  achievable  during  water-cooled  copper  mold 
casting  of  0.5  mm  and  1  mm,  respectively.2,4  The  two  other 
alloys,  Ca60Mg15Zn25  and  Ca60Mg20Zn20  are  better  bulk  glass 
formers  and  their  maximum  fully  amorphous  plate  thick- 


TABLE  I.  Density  (in  g/cm3  and  atoms /A3)  of  Ca-Mg-Zn 
amorphous  alloys  produced  in  this  work. 


Alloy 

Density 

g/cm3 

Atoms /A3 

Ca60Mg10Zn3o 

2.4481  ±0.0058 

0.031985 

Ca60Mg15Zn25 

2.2890  ±0.0073 

0.031301 

Ca60Mg20Zn20 

2.1499  ±0.0024 

0.030837 

Ca6oMg25Zn15 

2.0043  ±  0.0044 

0.030227 

nesses  were  identified  to  be  6  mm  and  4  mm,  respectively.  In 
order  to  exclude  the  effects  of  casting  conditions,  amorphous 
specimens  of  all  four  alloys  were  prepared  by  melt  spinning 
in  the  form  of  ~0.2-mm-thick  ribbons.  The  densities,  p0,  of 
the  amorphous  alloys  were  measured  with  a  helium  pycnom¬ 
eter  AccuPyc  1330  VI. 03  and  the  values  (in  g/cm3  and 
atoms  /A3)  are  given  in  Table  I.  The  following  formula  was 
used  for  the  density  conversion: 

Po  (atoms/ A3)  =  WA/±  £'A)p0  (g/cm3),  (1) 

where  NA  is  the  Avogadro’s  number  and  ct  and  At  are  the 
atomic  fraction  and  atomic  mass,  respectively,  of  the  element 
i  (Ca,  Mg,  or  Zn)  in  a  given  alloy. 

B.  Neutron  and  x-ray  diffraction 

Neutron-diffraction  experiments  were  conducted  at  room 
temperature  under  vacuum  using  the  general  materials 
(GEM)  diffractometer  at  the  ISIS  high-intensity  pulsed  neu¬ 
tron  source  (Rutherford  Appleton  Laboratory,  Didcot, 
U.K.)  23  GEM  has  eight  detector  banks  which  collects  data 
over  a  wide  range  in  Q  (from  0.1  to  100  A-1),  where  Q 
=  477  sin  0/X  is  the  magnitude  of  the  scattering  vector  for  a 
neutron  of  wavelength  X  scattered  at  an  angle  20.  The 
samples  were  loaded  in  10.3  mm  diameter  cylindrical  vana¬ 
dium  cans,  and  the  ribbons  were  crushed  to  allow  efficient 
packing  of  the  material.  The  thickness  of  the  can  walls  was 
25  jum  to  minimize  corrections.  The  exposure  time  for  each 
sample  was  about  3.5  h.  Due  to  the  high  degree  of  structural 
disorder  in  these  samples,  there  were  no  apparent  oscillations 
in  the  diffraction  pattern  beyond  25  A-1.  X-ray  diffraction 
experiments  were  conducted  at  room  temperature  on  a  Pana- 
lytical  Xpert-Pro  diffractometer  (ISIS  Facility)  using  Ag  Ka 
radiation.  The  samples  were  powdered  and  loaded  into  1  mm 
diameter  silica  capillaries  of  wall  thickness  0.01  mm,  which 
were  sealed  with  wax.  The  Q  range  extended  from  0.5  to 

18  A"1. 

ISIS  developed  software,  gudrun,24  gudrunx,25  and  AT¬ 
LAS  (Ref.  26)  were  used  to  reduce  and  correct  the  neutron 
and  x-ray  diffraction  data.  This  yields  the  experimental  total 
scattering  structure  factors,  SN(Q)  and  SX(Q),  and  the  total 
radial  distribution  functions  (RDFs)  GN(r)  for  the  four  stud¬ 
ied  Ca-Mg-Zn  metallic  glasses.  The  superscript  N  or  X  indi¬ 
cates  that  the  data  was  obtained  by  neutron  or  x-ray  diffrac¬ 
tion,  respectively.  These  functions  are  related  to  the  partial 
RDFs  (PRDFs),  g^(r),  by  the  following  equations:27 
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sn(Q)  ■■ 


—  2c,cmP 

Q  Q=x  1  JJo 


dgij(r)  ~  l]sin (Qr)dr,  (2) 


SX(Q)  =  ~^T  2  Cfj  f  r[g..(r )  _  l]sin (Qr)dr, 

.Jo 


Q  U— 1  ~  0 

2  MQ)2 


k=l 


(3) 


3 

GN(r )  =  2  CiCjbibjigijir)  -  1],  (4) 

0=1 

In  Eqs.  (2)-(4),  q,  bb  and  /)(g)  are,  respectively,  the  atomic 
fraction,  coherent  bound  neutron-scattering  length,  and  x-ray 
scattering  factor  of  element  i  (Ca,  Mg,  or  Zn),  p0  is  the 
average  number  density  of  the  material  (in  atoms /A3),  and  r 
is  a  distance  from  an  average  origin  atom  point  in  the  amor¬ 
phous  structure.  The  GN(r)  functions  for  the  four  studied 
Ca-Mg-Zn  glasses  were  obtained  by  Fourier  transformation 
of  the  respective  SN(Q)  functions  using  the  Lorch  modifica¬ 
tion  function28  with  a  maximum  momentum  transfer  <2max  of 

25  A-1. 


C.  Reverse  Monte  Carlo  simulation 

Three-dimensional  models  of  the  amorphous  structures  of 
the  Ca-Mg-Zn  alloys  that  agree  with  the  experimental  data, 
alloy  density  and  closest  approach  constraints,  were  gener¬ 
ated  with  the  aid  of  the  RMC  simulation  technique  described 
in  detail  in  Ref.  29.  Simulation  boxes  with  periodic  boundary 
conditions  contained  5400  atoms,  and  the  box  size  was  de¬ 
termined  by  the  density  of  the  material  being  modeled.  The 
input  data  used  in  each  simulation  were  the  experimental 
neutron  and  x-ray  total  structure  factors,  SN(Q)  and  SX(Q), 
and  the  radial  distribution  function,  GN(r).  Gx(r )  was  not 
used  because  a  noticeable  error  in  calculating  this  function 
from  SX(Q),  determined  in  a  narrow  Q  range,  was  expected. 
Indeed,  the  g-dependent  x-ray  scattering  factors, /)(<2),  make 
corrections  of  SX(Q)  and,  especially,  Gx(r )  from  the  experi¬ 
mental  x-ray  data  more  difficult,  and  so  the  results  are  less 
reliable,  than  the  neutron  data.30  Use  of  both  SN(Q)  and 
GN(r)  and  only  SX(Q)  thus  gave  more  weight  to  the  more 
reliable  neutron  data. 

The  closest  approach  constraints  (i.e.,  minimum  inter¬ 
atomic  distances)  together  with  the  fixed  alloy  density  have 
been  found  to  improve  the  separation  of  partials  in  cases 
where  the  separation  matrix  is  poorly  conditioned.29,31  These 
two  main  constraints  were  held  throughout  the  whole  simu¬ 
lation  runs  and  severely  limited  the  number  of  structures  that 
were  consistent  with  the  experimental  data.  The  minimum 
pair  bond  distances,  r*™n,  between  i  and  j  atoms  were  com¬ 
puted  from  the  results  of  deconvolution  of  the  first  peak  of 
the  total  correlation  function,  TN(r),  with  Gaussian  functions, 
using  a  peakfit  version  4.12  program  (SeaSolve  Software, 
Inc.,  San  Jose,  CA),32 


T%j  =  4tti'p0  2  CjCjbfbjgJr) 


ij=  1 


3  3 

i=l  j^i  \27T(T ij 


exp 


(5) 


Here,  r and  oq  are,  respectively,  the  most  probable  (mode) 
bond  distance  and  the  standard  deviation  of  the  i-j  partial 
(Gaussian)  peak,  atj  is  the  peak  area,  and  r™in  =  rij-3crij.  Thus 
the  following  minimum  pair-bond  distances  were  applied 
in  the  RMC  simulations:  r^Ca=3.26  A,  r£anM  =2.96  A, 

min  _o  /co  X  r1™1  -?  SR  A  rmin  -9  20  A  and  rmin 
rCa-Zn“Z-OZ  'Mg-Mg-^5  rMg-Zn-Z-JU  dnQ  rZn-Zn 

=  2.17  A.  (In  this  paper,  we  use  the  term  “bond  distance”  for 
convenience  when  discussing  two  atoms  which  are  nearest 
neighbors,  and  this  should  not  be  taken  as  implying  that  there 
is  any  particular  type  of  bonding  between  the  two  atoms.) 
During  the  initial  stages  of  the  RMC  simulation,  these 
Gaussian  pair-bond  distribution  profiles  were  used  as  addi¬ 
tional  constraints  in  order  to  prevent  formation  of  physically 
unrealistic  structures  with  sharp  peaks  of  some  g^(r)  at  r™m. 
As  soon  as  a  stationary  difference  between  the  calculated  and 
experimental  RDFs  was  achieved,  these  additional  con¬ 
straints  were  removed  and  the  simulated  structure  was  fur¬ 
ther  refined  to  achieve  the  best  fit  between  the  experimental 
and  calculated  RDFs  and  total  structure  factors.  More  than 
10  million  accepted  random  atom  moves  occurred  for  each 
sample  before  the  final  configurations  were  collected,  which 
provided  total  independence  of  the  modeled  structures  on  the 
initial  configurations  of  randomly  distributed  atoms. 

The  RMC  computed  amorphous  structures  were  analyzed 
with  the  use  of  freeware  programs33  and  programs  developed 
by  coauthors  from  John  Hopkins  University.  These  programs 
allowed  the  calculation  of  structural  features  such  as  PRDFs; 
type  and  distribution  of  coordination  polyhedra  (atomic  clus¬ 
ters);  neighbor  coordination  numbers  (CNs)  and  triplet  angle 
correlations. 


III.  RESULTS 
A.  Diffraction  data 

Figure  3  shows  the  GN(r)  data  for  the  four  samples.  Sys¬ 
tematic  changes  with  composition  can  be  observed.  In  par¬ 
ticular,  the  addition  of  Mg  at  the  expense  of  Zn  results  in  a 
decrease  in  intensity  in  the  region  of  2. 8 -3. 3  A,  and  an 
increased  intensity  between  3.3  and  4.0  A.  Using  the  metal¬ 
lic  and  covalent  atomic  radii34,35  (Table  II)  to  calculate  inter¬ 
atomic  distances,  it  is  predicted  that  the  former  region  is 
comprised  of  overlapping  correlations  from  Zn-Zn,  Zn-Mg, 
and  Zn-Ca  pairs  while  the  latter  arises  from  Mg-Mg,  Mg-Ca, 
and  Ca-Ca  correlations.  The  relative  change  in  intensity 
across  the  first  peak  manifold  can  therefore  be  understood  as 
relating  directly  to  the  exchange  of  magnesium  atoms  for 
zinc  in  the  glass  structure.  However,  information  on  the  pre¬ 
cise  interatomic  distances  and  coordination  numbers  of  each 
interatomic  pairing  cannot  be  extracted  from  the  data  due  to 
the  close  overlapping  of  the  six  individual  pair  correlations 
in  this  system.  Such  detailed  information  on  the  PRDFs  can 
only  be  extracted  by  careful  modeling  of  the  data. 
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FIG.  3.  (Color  online)  Experimental  and  RMC  simulated  neu¬ 
tron  total  radial  distribution  functions  GN(r)  for  Ca60MgxZn40_x 
amorphous  alloys. 


FIG.  4.  (Color  online)  Experimental  and  RMC  simulated  x-ray 
and  neutron  total  scattering  structure  factors  S(Q )  for  Ca60Mg2oZn2o 
and  Ca60Mg25Zn15  amorphous  alloys. 


B.  Pair  correlations 

The  eventual  matches  between  the  simulated  and  experi¬ 
mental  GN(r )  for  the  four  Ca-Mg-Zn  amorphous  alloys,  after 
^107  RMC  steps,  are  shown  in  Fig.  3  and  the  total  scattering 
factors  are  shown  for  two  of  the  samples  in  Fig.  4.  An  ex¬ 
cellent  fit  indicates  that  the  RMC- simulated  amorphous 
structure  describes  the  experimental  data  very  well.  The  first 
GN(r )  peak  located  in  the  r  range  between  2.21  and  4.67  A, 
corresponds  to  distances  between  atoms  within  the  first  co¬ 
ordination  shell.  The  six  PRDFs  for  the  simulated  amorphous 
structure  of  Ca6oMg20Zn20  are  shown  in  Fig.  5.  The  PRDFs 
for  the  Ca  correlations  are  well  determined  due  to  the  abun¬ 
dance  of  Ca  in  the  glass.  A  single  g^-(r)  peak  inside  the  first 
coordination  shell  range  is  clearly  seen  in  the  Ca-Ca,  Ca-Mg, 
and  Ca-Zn  PRDFs,  and  the  bond  distances,  rljp  between  the 
Ca-Ca,  Ca-Mg,  and  Ca-Zn  pairs  corresponding  to  the  maxi¬ 
mum  of  this  peak  are  identified  to  be  3.84  A,  3.48  A,  and 
3.21  A,  respectively.  On  the  other  hand,  the  Mg-Mg,  Mg- 
Zn,  and  Zn-Zn  correlations  are  less  well  defined  due  to  the 
low  abundance  of  Mg  and  Zn  atoms  in  this  glass  system.  The 
first  peaks  in  the  PRDFs  for  these  correlations  are  noisy  and 
rather  wide.  It  is  likely  that  a  spectrum  of  bond  distances  is 
typical  for  the  Mg-Mg,  Mg-Zn  and,  especially,  Zn-Zn  corre¬ 
lations  within  the  first  shell  (i.e.,  inside  the  first  RDF  peak). 
A  similar  behavior  of  the  PRDFs  has  also  been  observed  for 
the  three  other  Ca-Mg-Zn  amorphous  alloys  in  this  study. 
The  interatomic  bond  distances  calculated  from  the  position 
of  the  maximum  intensity  of  the  first  peak  of  the  respective 


TABLE  II.  Metallic,  rm,  (Ref.  34)  and  covalent,  rc,  (Ref.  35) 
bond  distances  (in  A)  between  Ca,  Mg,  and  Zn  atom  pairs. 


Ca-Ca 

Ca-Mg 

Ca-Zn 

Mg-Mg 

Mg-Zn 

Zn-Zn 

rm  (A) 

3.94 

3.57 

3.31 

3.20 

2.94 

2.68 

rc  (A) 

3.52 

3.17 

2.98 

2.82 

2.63 

2.44 

PRDFs  for  the  four  Ca-Mg-Zn  amorphous  alloys  are  given  in 
Table  III. 

Coordination  number  weighted  bond  distances,  r™el8hted, 
between  atom  pairs  within  the  first  peak  of  the  total  RDFs  of 
the  Ca-Mg-Zn  amorphous  alloys  are  given  in  Table  IV.  They 
are  calculated  using  the  following  equation: 


where  rx  and  r2  are  the  minimum  and  maximum  cutoffs  for 
the  first  shell.  When  the  values  in  Table  IV  are  compared 
with  respective  values  in  Table  III,  it  is  found  that  rf.ei8hted 
are  always  higher  than  the  rtj  values  at  the  first  peak  maxi¬ 
mum  intensity.  The  difference  increases  with  a  decrease  in 
the  size  of  pair  atoms.  This  finding  indicates  that  the  gij(r) 
peaks  are  asymmetrical  (non-Gaussian),  especially  for  the 
Mg-Mg,  Mg-Zn,  and  Zn-Zn  atom  pairs,  which  is  an  indica¬ 
tion  that  these  pairs  have  a  wide  spectrum  of  bond  distances 
within  the  first  shell. 


IV.  DISCUSSION 
A.  Pair  correlations 

The  interatomic  bond  distances  calculated  from  the  posi¬ 
tion  of  the  maximum  intensity  of  the  first  peak  of  the  respec¬ 
tive  PRDFs  for  the  four  Ca-Mg-Zn  amorphous  alloys  (Table 
III)  are  about  2.3 -4.1  %  smaller  than  the  respective  metallic 
bond  distances,  rm,  and  about  6.5 -9.8  %  larger  than  the  re¬ 
spective  covalent  bond  distances,  rc  (see  Table  II).  The 
Ca-Zn  bonds  are  noticeably  shorter  and  Zn-Zn  bonds  are 
longer  in  marginal  glasses  (Ca60Mg10Zn30  and 
Ca60Mg25Zn15)  than  in  bulk  glass  formers  (Ca60Mg15Zn25 
and  Ca60Mg20Zn20).  Table  V  provides  characteristic  first- 
shell  bond  distances  between  Ca,  Mg,  and  Zn  atom  pairs  in 
several  crystalline  binary  intermetallics.36  There  are  discrete 
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FIG.  5.  RMC-simulated  partial  radial  distribution  functions  gy(r)  for  Ca-Ca,  Ca-Mg,  Ca-Zn,  Mg-Mg,  Mg-Zn,  and  Zn-Zn  atomic  pairs  in 
Ca60Mg2oZn2o  amorphous  alloy. 


values  of  the  bond  distances  constrained  by  long-range  order 
of  a  crystal  structure  in  these  intermetallics.  The  minimum 
bond  distances  between  Ca-Ca,  Ca-Zn,  Mg-Mg,  and  Zn-Zn 
pairs  in  the  intermetallic  phases  are  smaller  while  those  be¬ 
tween  Ca-Mg  and  Zn-Mg  pairs  are  larger,  than  the  lengths  of 
respective  metallic  bonds  (compare  the  data  in  Tables  II  and 
V)  indicating  presence  of  covalent  components  in  these  crys¬ 
tal  structures.  The  peak  Ca-Mg  and  Zn-Mg  pair  distances  in 
the  amorphous  Ca-Mg-Zn  alloys  given  in  Table  III  are 
smaller  than  the  minimum  distances  between  respective  pairs 
in  the  crystalline  intermetallic  phases  while  the  peak  rtj  val¬ 
ues  for  other  four  pairs  in  the  amorphous  alloys  are  close  to 
the  minimum  r values  in  the  crystalline  intermetallics.  This 
comparative  analysis  indicates  that  the  absence  of  long-range 
order  constraints  allows  shortening  of  interatomic  distances 
between  the  nearest  neighbor  atoms  in  the  amorphous  struc¬ 
ture,  as  compared  to  the  binary  intermetallics  and  pure 
metals. 


The  observation  that  some  nearest-neighbor  bond  lengths 
are  shorter  than  those  in  the  competing  long-range  ordered 
crystals  may  help  explain  why  these  glasses  are  easy  glass 
formers.  The  equilibrium  crystalline  phases  that  compete 
with  the  amorphous  phase  appear  to  minimize  total  energy 
by  sacrificing  optimal  short-range  order  to  achieve  beneficial 
long-range  order.  The  absence  of  long-range  atomic  order  is 
an  energetic  disadvantage  but  allows  metallic  glasses  to 
maintain  optimized  short-range  atomic  interactions  that  help 
reduce  the  energy  difference  between  the  metastable  glass 
and  equilibrium  crystal(s).  The  largest  contribution  to  the 
stability  of  condensed  phases  is  expected  to  come  from 
nearest-neighbor  interactions  so  that  the  optimized  short- 
range  interactions  in  metallic  glasses  can  produce  structures 
that  have  only  a  small  energetic  disadvantage  relative  to  the 
equilibrium  structure.  Kinetic  constraints  from  quenching  re¬ 
strict  the  long-range  atomic  redistribution  needed  to  achieve 
long-range  order  and  to  further  minimize  the  total  system 


TABLE  III.  The  interatomic  bond  distances  (r^,  in  A),  that  correspond  to  the  position  of  the  maximum 
intensity  of  the  first  peak  in  the  respective  PRDFs  of  Ca60MgxZn40_x  amorphous  alloys.  The  bond  distance 
values  averaged  for  the  four  alloys,  (r^)aver,  as  well  as  the  deviations  of  these  values  from  metallic,  rm,  and 
covalent,  rc,  bond  distances  (see  Table  II),  are  also  given  here. 


rij 

(A) 

r  Ca-Ca 

r Ca-Mg 

r  Ca-Zn 

^Mg-Mg 

r Mg-Zn 

r Zn-Zn 

Ca60Mg10Zn30 

3.82 

3.44 

3.14 

3.02 

2.84 

2.66 

Ca60Mg|5Zn25 

3.80 

3.48 

3.19 

3.06 

2.78 

2.56 

Ca60Mg2oZn2o 

3.84 

3.48 

3.21 

3.02 

2.84 

2.56 

Ca60Mg25Zni5 

3.81 

3.52 

3.16 

3.09 

2.85 

2.63 

v  */)  aver 

3.82 

3.48 

3.18 

3.05 

2.87 

2.60 

I00%(raver/rm-l) 

-3.1 

-2.5 

-4.1 

-4.8 

-2.3 

-2.9 

100%(raver/rc-l) 

8.5 

9.8 

6.5 

8.1 

9.2 

6.7 
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TABLE  IV.  Weighted  bond  distances,  r^ei8hted  (in  A),  between  the  atom  pairs  within  the  first  global  peak 
of  the  corresponding  PRDFs  of  the  Ca60MgxZn40_x  amorphous  alloys  [Eq.  (6)].  The  weighted  bond  distance 
values  averaged  for  the  four  alloys,  (r^ei8hted)ayc[,  as  well  as  the  difference  between  (r^-eighted) ^ r  and  (r^)aver, 
(taken  from  Table  III)  are  also  given. 


weighted 

(A) 

Ca-Ca 

Ca-Mg 

Ca-Zn 

Mg-Mg 

Mg-Zn 

Zn-Zn 

Ca60Mg10Zn3o 

3.83 

3.55 

3.32 

3.35 

3.20 

3.21 

Ca60Mg  |5Zn25 

3.85 

3.57 

3.29 

3.40 

3.33 

3.30 

Ca60Mg2oZn2o 

3.85 

3.53 

3.35 

3.44 

3.34 

3.36 

Ca60Mg25Zn|5 

3.85 

3.52 

3.43 

3.43 

3.19 

3.38 

(  yweighted\ 

\'ij  /aver 

3.84 

3.54 

3.35 

3.40 

3.26 

3.31 

(  weight ed\  /  \ 

v  ij  )  aver  V  ij)  aver 

0.02 

0.06 

0.17 

0.35 

0.39 

0.71 

energy,  thus  favoring  glass  formation.  The  present  results 
suggest  that  metallic  glasses  not  only  have  SRO  but  may 
have  “better”  short-range  topological  order  (in  terms  of  op¬ 
timal  bond  length  and/or  atomic  arrangement)  than  the  com¬ 
peting  crystals. 

The  Ca-Mg-Zn  amorphous  structures  show  asymmetrical 
(non-Gaussian)  distribution  of  the  neighbor  bond  lengths, 
which  result  in  weighted  (i.e.,  mean)  bond  distances, 
rweighted ,  between  at0m  pairs  to  be  always  higher  than  the 
most  probable  (i.e.,  mode)  pair  distances,  rtj,  that  correspond 
to  the  position  of  the  first  peak  maximum  intensity  in  the 
respective  PRDFs  (see  Tables  III  and  IV).  The  difference 
between  rJ.eishted  anc[  r..  [s  minimum  for  Ca-Ca  pairs 

TABLE  V.  Characteristic  bond  distances  for  coordination  poly- 
hedra  in  several  crystalline  binary  intermetallic  phases  (Ref.  36). 


r 


(A)  rCa-Ca  rCa-Mg  rCa-Zn  rMg-Mg  r Mg-Zn  r Zn-Zn 


CaMg2 

3.81 

3.62 

3.05 

3.13 

3.18 

Ca3Zn 

3.78 

3.18 

4.15 

3.88 

3.56 

3.91 

4.17 

3.95 

4.15 

4.32 

CaZn 

3.90 

3.18 

2.62 

4.05 

3.22 

4.20 

MgZn2 

3.17 

3.04 

2.53 

2.61 

2.64 

Ca5Zn3 

3.61 

3.14 

2.69 

3.76 

3.29 

4.07 

3.34 

4.19 

Minimum 

3.61 

3.62 

3.14 

3.05 

3.04 

2.53 

Maximum 

4.32 

3.62 

4.17 

3.18 

3.04 

4.15 

(—0.02  A),  increases  in  the  order  of  Ca-Mg,  Ca-Zn,  Mg- 
Mg,  Mg-Zn,  and  is  maximum  for  Zn-Zn  pairs 
(—0.71  A).  The  clearly  identified  and  more  symmetrical 
PRDF  peaks  for  the  Ca-Ca,  Ca-Mg,  and  Ca-Zn  atom  pairs 
within  the  first  coordination  shell,  with  the  mean  and  mode 
interatomic  distances  equal  to  or  shorter  than  the  neighbor 
bond  distances  found  in  crystalline  metals,  may  suggest  that 
the  atoms  in  these  pairs  attract  and  are  touching  each  other. 
On  the  other  hand,  the  highly  asymmetrical  distributions  of 
the  bond  distances  for  the  Mg-Mg,  Mg-Zn,  and  Zn-Zn  atom 
pairs  is  an  indication  that  these  pairs  have  a  spectrum  of 
bond  distances  within  the  first  shell.  Although  some  of  these 
pairs  are  close  packed  and  their  mode  pair  distances  are 
shorter,  their  mean  distances  are  much  longer  than  the  re¬ 
spective  metallic  bond  distances  (compare  Tables  II-IV),  in¬ 
dicating  that  Mg  and  Zn  are  less  attractive  to  each  other  than 
to  Ca.  Based  on  these  results,  one  can  suggest  that,  in  aver¬ 
age,  Mg-  and  Zn-centered  clusters  prefer  to  have  Ca  atoms  as 
nearest  neighbors,  while  Mg  and  Zn  atoms  in  the  first  shell 
just  fill  the  gaps  between  the  Ca  atoms,  so  that  their  distances 
from  the  center  atom  (Mg  or  Zn)  vary  depending  on  the  gap 
sizes.  These  observations  also  suggest  that  pairs  of  Zn  atoms 
more  strongly  prefer  to  be  separated  than  other  atom  pairs. 

B.  Coordination  polyhedra 

Atoms  in  metallic  glasses  tend  to  arrange  themselves  so 
as  to  maximize  the  local  packing  density20,37-39  and  the  finite 
sizes  of  atoms  impose  geometrical  restrictions  on  the  pos¬ 
sible  local  atom  environments.40,41  The  local  environment  of 
an  atom  is  generally  treated  in  terms  of  the  geometry  of  the 
shell  formed  by  its  nearest- neighbor  atoms  42,43  The  chemis¬ 
try  (i.e.,  the  number  of  different  atom  species  and  their  rela¬ 
tive  positions)  of  the  first  coordination  shell  provides  addi¬ 
tional  important  information  about  the  local  environment. 
Bonding  in  metallic  glasses  is  generally  nondirectional,  and 
the  nearest  neighbors  for  any  chosen  atom  i  are  often  defined 
as  those  atoms  that  are  either  in  contact  with  /,  or  separated 
from  it  by  a  distance  smaller  than  a  certain  threshold  value. 
The  latter  is  defined  globally  as  the  position  of  the  first  mini¬ 
mum  of  the  RDF  of  the  whole  atomic  array.  Therefore,  this 
approach  is  not  sensitive  to  the  details  of  the  first  shell  ge¬ 
ometry  and/or  to  local  density  and  chemistry  fluctuations. 
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A  more  advanced  approach  is  to  define  the  nearest  neigh¬ 
bors  as  atoms  that  have  common  faces  in  their  Voronoi 
polyhedra.37  This  definition  of  the  first  coordination  shell  at¬ 
oms  is  straightforward  and  reflects  local  details  of  the  atomic 
packing  topology.  In  this  approach,  a  coordination  polyhe¬ 
dron  with  vertices  located  in  positions  of  the  first  coordina¬ 
tion  shell  atoms  and  with  edges  coinciding  with  the  inter¬ 
atomic  bonds  is  defined  for  any  chosen  atom  / 43  Each 
coordination  polyhedron,  which  is  also  called  an  /-centered 
cluster,15,20  is  associated  with  the  respective  /-centered 
Voronoi  polyhedron  and  can  therefore  be  assigned  the  char¬ 
acteristic  Voronoi  signature  (v3,v4,v5,v6).  For  the  Voronoi 
polyhedron,  vm  is  the  number  of  faces  containing  m  edges 
while  for  the  respective  coordination  polyhedron,  vm  is  the 
number  of  the  vertices  common  to  m  polyhedron  edges  (or 
faces).44  In  the  latter  case,  m  is  also  called  the  vertex  coor¬ 
dination.  The  bond  connecting  the  center  atom  /  with  an 
m-coordinated  neighbor  atom/vertex  j  is  perpendicular  to  an 
m-edged  face  in  the  respective  /-centered  Voronoi  polyhe¬ 
dron,  and  thus  called  an  ra-fold  bond.  Those  /-centered  clus¬ 
ters  with  the  same  Voronoi  signature  are  considered  to  be 
topologically  equivalent  (even  though  they  may  not  be  iden¬ 
tical)  because  they  can  be  transformed  into  each  other  by 
‘elastic’  deformation  without  changing  the  number  of  verti¬ 
ces  and  connecting  edges.  In  addition  to  assigning  topology, 
the  Voronoi  signature  also  defines  the  total  CN  of  the  cluster 
as  CN=Sum.  Clusters  with  the  same  Voronoi  signature  can 
however  be  chemically  different  if  they  consist  of  different 
elements.  Therefore,  in  addition  to  the  Voronoi  signature, 
partial  coordination  numbers,  i.e.,  the  number  of  atoms  of 
different  species,  should  also  be  known  for  a  more  complete 
description  of  the  SRO  in  the  amorphous  structure. 

Figure  6  illustrates  the  types  and  fractions  of  Ca-,  Mg- 
and  Zn-centered  coordination  polyhedra  found  in  the  RMC- 
simulated  amorphous  structures  of  the  four  Ca-Mg-Zn  al¬ 
loys.  Although  many  types  of  coordination  polyhedra  are 
present  in  the  amorphous  structures,  the  most  common  clus¬ 
ters  can  be  identified.  These  are  (0,1,10,2)  for  Ca-centered, 
(0,2,8, 1)  for  Mg-centered,  and  (0,3, 6,0),  (0,2,8, 1),  (0,3,6, 1), 
and  (0,2, 8,0)  for  Zn-centered  clusters.  Three  examples  of  the 
clusters  extracted  from  the  amorphous  structures  are  illus¬ 
trated  in  Fig.  7.  The  characteristic  coordination  polyhedra  in 
all  these  Ca-Mg-Zn  glasses  are  not  dominated  by  (0,0,12,0) 
icosahedra  and  Mg-  and  Zn-centered  polyhedra  with  CNs  of 
11,  10,  and  9  are  actually  more  common.  This  is  similar  to 
the  Cu-Zr  BMGs  on  the  Zr-rich  side,  where  solute  Cu  atoms 
are  surrounded  mostly  by  solvent  Zr  and  have  CNs  lower 
than  12,  but  different  from  the  more  Cu-rich  Cu-Zr  (and 
Cu-Zr-Al)  BMGs,  where  the  Cu-centered  (0,0,12,0)  icosahe¬ 
dra  with  both  Cu  and  Zr  in  the  shell  are  the  dominant  struc¬ 
tural  unit  influencing  the  GFA  and  properties.11-13 

Whereas  the  simulated  amorphous  structures  of 
Ca60MgxZn40_x  alloys  contain  many  types  of  coordination 
polyhedra,  the  competing  crystal  phases,  CaMg2  and  CaZn,10 
contain  only  three  and  two  types  of  polyhedra,  respectively.36 
Figure  8  shows  the  three  principal  clusters  (one  is  Ca- 
centered  and  two  are  Mg-centered)  that  are  present  in  the 
CaMg2  intermetallic  phase,  and  Fig.  9  shows  the  two  clusters 
(Ca-  and  Zn-centered,  respectively)  found  in  the  CaZn  inter¬ 
metallic  phase.  The  Ca-centered  clusters,  (0,0,12,4)  in 
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FIG.  6.  (Color  online)  Types  of  Ca-,  Mg-,  and  Zn-centered  clus¬ 
ters  and  their  fractions  in  four  Ca60MgxZn40_z  amorphous  alloys. 

CaMg2  and  (0,1,10,6)  in  CaZn,  are  not  present  in  the  ana¬ 
lyzed  amorphous  structures  where  the  number  of  atoms 
around  each  Ca  atom  is  found  to  be  significantly  less  than 
16-17,  indicating  that  the  local  environment  of  Ca  is  not 
efficiently  packed  in  the  amorphous  structure.  Both  Mg- 
centered  clusters  in  CaMg2  are  of  the  (0,0,12,0)  type  and  the 
Zn-centered  cluster  in  CaZn  is  (0,3, 6,0).  Similar  Mg-  and 
Zn-centered  clusters  are  also  common  in  the  amorphous 
structures  (see  Fig.  6).  This  is  different  from  transition  metal 
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FIG.  7.  (Color  online)  Configurations  of  three  typical  (a)  Ca- 
centered  (0,1,10,2),  (b)  Mg-centered  (0,2,8, 1),  and  (c)  Zn-centered 
(0,3,6, 1)  clusters  in  the  Ca60MgxZn40_x  amorphous  alloys. 

Zr-Cu  glasses,  where  the  dominant  icosahedra  and  related 
clusters  are  not  found  in  the  competing  crystalline  com¬ 
pounds. 

Of  particular  note  is  the  increase  in  the  fraction  of 
(0,0,12,0)  Mg  clusters  from  —0.03  to  —0.13  when  the  Mg 
content  is  increased  from  20  to  25  at.  %.  This  increase  re¬ 
sults  in  the  cluster  being  the  third  most  common  in  the  glass, 


FIG.  8.  (Color  online)  (a)  Ca-centered  (0,0,12,4)  and  (b)  and  (c) 
Mg-centered,  both  are  (0,0,12,0),  coordination  polyhedra  in  the 
CaMg2  crystal  structure  (Ref.  36). 


FIG.  9.  (Color  online)  (a)  Ca-centered  (0,1,10,6)  and  (b)  Zn- 
centered  (0,3, 6,0)  coordination  polyhedra  in  the  CaZn  crystal  struc¬ 
ture  (Ref.  36). 

with  the  fraction  of  the  two  most  common  clusters,  (0,2,8, 1) 
and  (0,2, 8, 2),  also  increasing  significantly  over  this  compo¬ 
sitional  range.  Furthermore,  in  contrast  to  the  other  clusters 
found  in  the  models,  these  latter  two  clusters  are  11  or  12 
coordinated  and  have  more  than  two  thirds  of  the  vertices 
with  the  coordination  number  m  =  5,  making  them  most  like 
the  (0,0,12,0)  icosahedral  cluster.  This  tendency  toward  12- 
coordinated  Mg  and  m  =  5  coincides  with  a  reduction  in  the 
maximum  amorphous  thickness  from  4  to  1  mm,  and  the 
formation  of  CaMg2.  While  less  obvious,  there  is  a  similar 
increase  in  the  number  of  Ca  atoms  with  coordination  num¬ 
bers  of  12  and  13  and  a  predominance  of  u5,  with  the  three 
most  dominant  clusters  for  the  Ca60Mg25Zn15  glass  being 
(0,1,10,2),  (0,3, 6, 4)  and  (0,2, 8, 2),  all  of  which,  in  addition  to 
(0,0,12,0),  are  present  in  significantly  higher  amounts  than  in 
the  Ca60Mg20Zn30  glass.  No  such  trends  toward  CaZn-like 
clusters  can  be  seen  in  the  Zn  environments  but  the  fraction 
of  “favorable”  (0,3, 6,0)  clusters  is  found  to  reduce  dramati¬ 
cally  as  the  Mg  content  increases  from  20  to  25  at.  %,  and 
CaMg2  is  formed  in  preference  to  CaZn.  Therefore,  increases 
in  the  fraction  of  local  clusters  that  are  topologically  similar 
to  crystalline  compounds  can  be  correlated  with  the  reduc¬ 
tion  in  GFA  but  their  presence  does  not  negate  glass  forma¬ 
tion.  This  is  true  for  oxide  glasses  such  as  amorphous  Si02 
and  for  the  binary  Pd-Si  BMG,45  and  the  BMG-forming  Ca- 
Mg-Zn  case  shows  that  it  is  also  true  for  more  complex  me¬ 
tallic  glasses. 

C.  Coordination  numbers 

The  distributions  of  the  atomic  clusters  by  the  total  CN  in 
the  first  shell  are  shown  in  Fig.  10.  There  is  no  notable  dif¬ 
ference  in  the  distributions  between  good  (X=  15  and  20)  and 
marginal  (X=10  and  25)  glass  formers.  The  most  common 
coordination  numbers  are  CN=13  for  the  Ca-centered  clus¬ 
ters,  CN=11  and  12  for  Mg-centered  clusters  and  CN=10 
for  Zn-centered  clusters. 

Figure  11  shows  the  average  number  of  Ca,  Mg,  and  Zn 
atoms  in  the  first  shell  of  the  Ca-,  Mg-,  and  Zn-centered 
clusters  in  the  Ca-Mg-Zn  amorphous  alloys.  Ca  atoms  are 
most  common  in  the  first  shell  of  the  clusters  and  an  increase 
in  the  concentration  of  Mg  from  10%  to  25%  almost  does  not 
affect  the  number  of  Ca  atoms  in  all  three  types  of  clusters. 
The  respective  Ca-Ca,  Mg-Ca,  and  Zn-Ca  partial  coordina¬ 
tion  numbers  (PCNs)  are  ^Ca_Ca=8.1-8.5,  ftMg_Ca=7.2-7.9, 
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FIG.  10.  (Color  online)  Distributions  of  Ca-centered,  Mg- 
cantered,  and  Zn-centered  clusters  by  their  total  coordination  num¬ 
bers  in  four  Ca60MgxZn40_z  amorphous  alloys. 


and  nZn_Ca=6J -1 .2.  In  comparison,  an  increase  in  the  con¬ 
centration  of  Mg  does  result  in  an  increase  in  the  number  of 
Mg  atom  neighbors  (i.e.,  ftCa_Mg  increases  from  1.2  to  3.0  and 
^Mg-Mg  increases  from  0.9  to  3.7)  and  a  corresponding  de¬ 
crease  in  the  number  of  Zn  atom  neighbors  (^ca-Zn  decreases 
from  3.3  to  1.8  and  nMg_Zn  decreases  from  2.3  to  0.6)  in  the 
Ca-  and  Mg-centered  clusters.  The  number  of  Zn  and  Mg 
atoms  in  the  first  shell  of  Zn-centered  clusters  does  not  vary 
linearly  and  a  link  between  the  GFA  and  the  local  environ¬ 
ments  can  be  made.  Namely,  marginal  glasses  (with  10%  and 
25%  Mg)  have  more  Zn  and  less  Mg  in  the  first  shell  of  the 
Zn-centered  clusters  than  good  glass-forming  alloys  (with 
15%  and  20%  Mg). 

The  CN  values  for  the  Zn-  and  Mg-centered  clusters  are 
higher  than  the  CN  values  of  9  and  10,  expected  for  these 
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FIG.  11.  (Color  online)  Partial  coordination  numbers  in 
Ca60MgxZn40_z  amorphous  alloys  as  a  function  of  Mg  content. 


clusters  from  the  solute-centered  ECP  model.20  However,  the 
ECP  model  assumes  that  only  solvent  (Ca)  atoms  are  present 
in  the  first  coordination  shell.  If  some  Ca  atoms  are  replaced 
by  the  smaller  Mg  and  Zn  atoms,  a  larger  number  of  atoms 
can  be  expected  in  the  first  shell.  This  idea  is  explored  quan¬ 
titatively  by  calculating  the  packing  efficiency  around  Ca, 
Mg,  and  Zn  atoms.  Using  the  partial  coordination  numbers  in 
Fig.  11,  atomic  radii  in  Table  III  and  the  analytic  approach  in 
Ref.  38  the  following  results  are  found.  Packing  around  Mg 
atoms  ranges  from  96%  for  Ca60Mg10Zn30  to  103%  for 
Ca60Mg15Zn25  and  Ca60Mg20Zn2o.  Within  an  anticipated  ac¬ 
curacy  of  ±10%,  these  results  show  that  all  of  the  available 
space  in  the  first  coordination  shell  of  Mg  atoms  is  fully 
occupied  by  the  combination  of  Ca,  Mg,  and  Zn  atoms  indi¬ 
cated  in  Fig.  11.  The  packing  efficiency  around  Zn  atoms  is 
similarly  efficient  and  ranges  from  103%  for  Ca60Mg10Zn30 
to  109%  for  Ca60Mg20Zn20  and  Ca60Mg25Zn15.  However,  the 
packing  around  Ca  atoms  ranges  from  87-88  %  and  does 
not  reaches  maximum  packing  efficiency.  This  result  for  Ca- 
Mg-Zn  glasses  is  similar  to  the  results  for  an  Al87Y8Ni5 
glass,38  where  packing  around  Y  and  Ni  solutes  was  within 
100  ±  10%  and  packing  around  the  solvent  (Al)  atoms  was 
only  69%.  A  notable  difference  is  that  packing  around  Ca 
solvent  atoms  in  the  Ca-Mg-Zn  BMGs  is  significantly  higher 
than  packing  around  the  Al  solvents  in  the  marginal 
Al87Y8Ni5  glass,  suggesting  that  the  more  efficient  packing 
of  solvent  atoms  around  Ca  contributes  to  the  higher  GFA 
relative  to  the  Al-Y-Ni  glass. 

A  local  packing  efficiency  that  exceeds  100%  around  Mg- 
and  Zn-centered  clusters  is  consistent  with  the  broad  range  of 
Zn-Zn,  Zn-Mg,  and  Mg-Mg  interatomic  separations  in  the 
first  shell.  Specifically,  the  packing  efficiency  reported  in  the 
preceding  paragraph  was  estimated  under  the  assumption 
that  the  first  shell  atoms  are  in  contact  with  the  central  atom. 
To  reduce  compression  strains  associated  with  “overpopula¬ 
tion”  in  the  first  coordination  shell,  some  Mg  and  Zn  atoms 
occupy  positions  farther  from  the  central  atom.  This  gives  a 
wider  and  less  symmetrical  first  PRDF  peak  for  Mg-Mg, 
Mg-Zn,  and  Zn-Zn  pairs.  As  a  result,  the  differences  between 
the  pair  distances  corresponding  to  the  maximum  peak  inten¬ 
sity  (mode)  and  weighted  average  (median)  become  large 
(see  Tables  III  and  IV).  On  the  other  hand,  for  Ca-centered 
clusters,  where  the  local  packing  efficiency  was  found  to  be 
—  87-88  %,  all  the  first  shell  atoms  contact  the  surface  of  the 
central  Ca  atom.  As  a  result,  a  very  sharp  first  peak  on  the 
corresponding  Ca-Me  PRDFs  (Me  represents  Ca,  Mg,  or  Zn) 
is  observed.  Moreover,  the  peak  and  weighted  distances  for 
the  Ca-Me  pairs  have  almost  the  same  values. 

As  it  has  been  mentioned  earlier,  the  neighbor  environ¬ 
ment  of  a  specific  j  atom  in  the  first  coordination  shell  of  an 
/-centered  cluster  can  be  described  by  the  vertex  coordina¬ 
tion  number  trip  i.e.,  the  number  of  neighbors  of  the  j  atom 
located  in  the  first  shell  of  the  /-centered  cluster39,43  or  the 
number  of  edges  (bonds)  entering  the  vertex  j.  Figure  12 
shows  the  fraction  of  3-,  4-,  5-  and  6-  coordinated  vertices  in 
the  Ca-,  Mg-,  and  Zn-centered  clusters.  In  all  three  cluster 
types,  the  vertex  coordination  m= 5  prevails.  The  dominance 
of  polyhedra  with  4-,  5-,  and  6-bond  vertices  has  been 
shown42,43  to  be  an  indication  of  dense  atomic  packing  of  the 
amorphous  structure.  Moreover,  a  high  fraction  of  5 -bond 
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FIG.  12.  (Color  online)  Fraction  of  vertices  with  3,  4,  5,  and  6 
neighbor  atoms  (bonds)  in  Ca-,  Mg-,  and  Zn-centered  clusters  in 
four  Ca60MgxZn40_x  amorphous  alloys. 


vertices  can  indicate  dominance  of  pentagonal  bipyramids  as 
structural  units  in  the  amorphous  structure.  Two  pentagonal 
bipyramids  with  a  common  fivefold  vertex  connected  by  ten 
tetrahedra  with  common  faces  form  a  coordination  icosahe¬ 
dron,  which  is  believed  to  be  the  most  common  cluster  in 
some  important  BMGs,  such  as  Cu-rich  Cu-Zr.46  In  Ca-Mg- 
Zn,  however,  many  of  the  pentagonal  bipyramids  do  not 
form  complete  (0,0,12,0)  icosahedra.  This  is  clearly  indi¬ 
cated  by  inspection  of  the  Voronoi  signatures,  (v3,v4,v5,v6), 
in  Fig.  6,  which  shows  that  none  of  the  values  of  u3,  u4,  or  v6 
are  larger  than  v5,  regardless  of  the  value  of  the  CN.  This  is 
also  indicated  by  the  low  fraction  of  clusters  with  the 
(0,0,12,0)  Voronoi  index. 

D.  Bond-angle  distributions 

Bond  angle  distributions  can  provide  additional  informa¬ 
tion  about  the  most  common  local  structures  around  cluster 
vertices  because  characteristic  vertex  angles  are  very  sensi¬ 
tive  to  the  vertex  type  (e.g.,  number  and  length  distributions 
of  vertex  bonds).  Figure  13  shows  Ca-Ca-Ca,  Ca-Mg-Ca, 
and  Ca-Zn-Ca  bond  angle  distributions  in  the  Ca-Mg-Zn 
amorphous  structures.  The  upper  limit  of  the  bond  lengths 
has  been  set  to  4.6  A,  which  corresponds  to  the  minimum 
after  the  first  peak  of  the  total  RDFs.  A  characteristic  feature 
for  all  four  alloys  is  the  first  peak  located  near  60°.  The  exact 
60°  angle  position  for  the  Ca-Ca-Ca  triplets  corresponds  to 
the  close  packing  of  three  equal  hard  spheres,  forming  equi¬ 
lateral  triangles  [Fig.  13(d)].  For  the  Mg-  and  Zn-centered 
triplets  (i.e.,  Ca-Mg-Ca  and  Ca-Zn-Ca),  the  first  peak  is  lo¬ 
cated  at  —64°  and  —68°,  respectively.  These  shifts  of  the 
first  peak  angle  to  higher  values  can  be  explained  by  the 
smaller  atomic  radii  of  Mg  and  Zn  atoms.  From  the  known 
bond  angle  values,  the  most  common  bond  distance  ratios 
are  estimated  to  be  rCa_Mg/rCa_Ca=0.94  and  rCa_Zn/rCa_Ca 
=  0.89,  which  leads  to  the  effective  atomic  radius  ratios  of 
(rMg/rca)eff=0.88  and  (/zn/rca)eff-0.78.  These  calculated  ra- 


(a)  ©  (degrees) 


(b)  ©(degrees) 


FIG.  13.  (Color  online)  (a)  Ca-Ca-Ca,  (b)  Ca-Mg-Ca,  and  (c) 
Ca-Zn-Ca  bond-angle  distributions  in  Ca60MgzZn40_x  amorphous 
alloys.  A  pentagonal  pyramid  with  two  characteristic  peak  bond 
angles,  a  —  60° -67°  and  /3~  107° -120°,  is  shown  in  (d)  and  rep¬ 
resents  the  most  common  configuration  in  the  first  coordination 
shell.  Me  in  figure  (d)  represents  Ca,  Mg,  or  Zn. 

dius  ratios  are  higher  than  both  the  metallic,  (rMg/rCa)m 
=  0.81  and  (rZn/rCa)m= 0.68,  and  the  covalent  radius  ratios, 
(rMg/rca)c=0.80  and  (/zn/rca)c=0.69,  calculated  from  Table 
II.  The  observed  Ca-Mg-Ca  and  Ca-Zn-Ca  triplet  configura¬ 
tions  can  be  explained  by  larger  shortening  of  Ca-Ca  bonds 
than  Ca-Mg  and  Ca-Zn  bonds  (relative  to  the  respective  me¬ 
tallic  bond  distances),  which  is  in  agreement  with  the 
weighted  bond  distance  data  reported  in  Table  IV. 
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TABLE  VI.  Structural  parameters  of  Ca60Mg15Zn25  obtained 
from  ab  initio  MD  simulation. 
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The  second  bond-angle  distribution  peak,  which  is  located 
near  107°  for  Ca-Ca-Ca  triplets,  corresponds  to  the  interior 
angle  of  108°  of  a  regular  pentagon  [see  Fig.  13(d),  angle  0]. 
A  slightly  smaller  value  is  probably  due  to  slight  distortion 
of  the  pentagon,  when  one  of  the  five  Ca  atoms  in  the  pen¬ 
tagon  vertices  is  replaced  by  Mg  and/or  Zn.  Due  to  the 
shorter  Ca-Mg  (or  Ca-Zn)  bonds,  the  interior  angle  at  the  Mg 
(Zn)  vertex  of  the  pentagon  increases  at  the  expense  of  other 
interior  angles,  the  total  sum  of  which  is  540°  in  the  penta¬ 
gon.  Indeed,  the  second  peak  for  Ca-Mg-Ca  and  Ca-Zn-Ca 
triplets  is  located  at  —116°  and  —120°,  respectively.  Having 
a  pentagon  consisting  of  four  Ca  atoms  and  one  Mg/Zn  atom 
and  using  these  values  of  the  interior  angle  at  the  Mg/Zn 
vertex,  the  interior  angles  at  the  other  four  Ca  vertices  are 
calculated  to  be  106°/ 105°.  To  obtain  the  average  angle  of 
107°  for  the  Ca-Ca-Ca  triplets,  it  should  be  assumed  that  one 
Mg/Zn  atom  replaces  every  10th/15th  Ca  atom  in  the  penta¬ 
gon  configurations  present  in  the  amorphous  structure. 

The  two  characteristic  peaks  in  the  bond-angle  distribu¬ 
tions  correspond  to  the  angles  at  the  vertices  of  equilateral 
triangles  and  pentagons  and  support  our  earlier  conclusion 
that  tetrahedra  and  pentagonal  bipyramids  are  the  most  com¬ 
mon  atomic  configurations  in  the  Ca60MgxZn40_x  amorphous 
alloys. 

E.  Ab  initio  molecular-dynamics  simulation 

To  investigate  whether  the  structure  obtained  using  RMC 
is  physically  stable  and  favored,  we  have  performed  ab  initio 
molecular-dynamics  (MD)  simulation  for  the  representative 
composition  Ca60Mg15Zn25.  The  cubic  supercell  with  200  at¬ 
oms  and  periodic  boundary  condition  was  equilibrated  at 
1000  K  for  100  ps,  and  quenched  to  300  K  at  —10  K/ps, 
using  VASP.47The  box  size  was  chosen  to  correspond  to  the 
experimental  density.  The  final  configuration  at  300  K  was 
analyzed  using  the  same  program  and  parameters,  and  the 
main  structural  features  are  in  agreement  with  the  RMC  re¬ 
sults;  see  Table  VI  for  the  bond  length  and  the  CNs.  The 
fraction  of  fivefold  bonds  is  also  dominant  (—60%). 

F.  Structural  description  via  efficient  cluster  packing 

Structural  descriptions  using  the  ECP  model  have  recently 
been  established  for  binary  metallic  glasses20  but  extension 
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to  ternary  glasses  has  not  yet  been  accomplished.  As  a  result, 
only  general  consistency  between  the  ECP  model  and  the 
present  results  are  possible.  The  most  significant  finding  of 
the  present  research  that  validates  the  underlying  physical 
principle  of  the  ECP  model  is  the  formation  of  a  structural 
scaffold  consisting  of  efficiently  packed  solute-centered  clus¬ 
ters  surrounded  primarily  by  Ca  atoms.  Earlier  concepts  sug¬ 
gest  that  the  largest  solute  dominates  the  structure  by  form¬ 
ing  the  solute-centered  clusters  which  form  the  structural 
scaffold.18,20  The  relative  size  between  the  a  solute  and  the 
solvent  defines  the  number  of  solvent  sites  that  surround 
each  a  site,  and  smaller  solutes  fit  in  [3  and  y  sites  situated  at 
the  center  of  octahedra  and  tetrahedra  of  the  primary  clusters, 
respectively.  Although  the  largest  solute  has  been  assumed  to 
occupy  a  sites  based  on  its  topological  potency  (it  bonds  to 
more  solvent  atoms  than  smaller  solutes),  it  is  also  possible 
that  smaller  solutes  may  first  form  clusters  with  the  solvent 
atoms.  For  example,  solutes  that  bond  more  strongly  to  the 
solvent  may  be  a  more  appropriate  designation  of  the  a  sol¬ 
ute.  Since  interatomic  bond  enthalpies  are  generally  not 
known  for  condensed  solids,  this  possibility  has  not  yet  been 
evaluated  in  detail.  The  relative  interaction  energies  between 
Ca,  Mg,  and  Zn  atoms  have  been  proposed  based  on  esti¬ 
mated  heats  of  mixing,48  and  show  that  the  Ca-Zn  interaction 
(-22  kJ/mol)  is  much  stronger  than  the  Ca-Mg  and  Mg-Zn 
interactions  (-6  kJ/mol  and  -4  kJ/mol,  respectively). 
Heats  of  mixing  give  a  global  representation  of  interatomic 
bond  enthalpies  that  includes  both  the  number  of  bonds 
formed  and  the  energies  of  bonds  between  two  atoms.  Both 
the  numbers  of  bonds  and  the  energies  of  bonds  are  likely  to 
be  composition  dependent  so  that  the  relative  bond  enthalp¬ 
ies  obtained  from  heats  of  mixing  should  be  used  with  cau¬ 
tion.  Nevertheless,  the  difference  between  the  Ca-Zn  interac¬ 
tion  and  the  Ca-Mg  and  Mg-Zn  interactions  is  sufficiently 
large  to  suggest  that  it  may  be  appropriate  to  consider  Zn  as 
the  a  solute.  The  stronger  preference  for  Zn  atoms  to  avoid 
contact  with  other  Zn  atoms  supports  this  assignment. 

In  the  present  work,  Mg  and  Zn  are  each  considered  sepa¬ 
rately  as  a  solutes.  In  the  former  case,  a  total  solute  atom 
fraction  of  0.261  is  required  to  fill  the  a ,  /3,  and  y  sites,  and 
in  the  latter  case,  a  total  solute  fraction  of  0.282  is  needed. 
The  ECP  model  predicts  that  no  solute- solute  nearest- 
neighbor  pairs  will  be  present  in  the  amorphous  structure 
with  solute  concentrations  below  these  levels,  and  that  the 
number  of  solute- solute  neighbor  pairs  will  increase  with 
increasing  solute  concentration  beyond  these  critical  values. 
Considering  either  Zn  or  Mg  as  the  a  solute,  the  total  solute 
content  in  the  investigated  alloys  is  sufficient  to  fill  all  the 
solute  sites,  with  the  remaining  fraction  occupying  solvent 
sites.  All  four  alloys  show  solute-solute  bonding  that  vary  in 
proportion  to  the  solute  concentration  (Fig.  11),  and  this 
gives  general  agreement  with  the  ECP  model.  More  convinc¬ 
ing  support  requires  quantitative  comparison  of  predicted 
PCNs  from  the  ECP  model  with  the  values  determined  here. 
This  requires  specific  assignment  of  the  Mg  and  Zn  solute 
atoms  to  a ,  /3,  and  y  solute  sites  and  to  solvent  sites.  Topo¬ 
logical  models  required  for  this  comparison  are  not  yet  avail¬ 
able  but  the  data  in  the  present  work  are  expected  to  give 
important  insights  into  development  of  such  models. 
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V.  CONCLUSIONS 

(1)  A  combination  of  neutron  and  x-ray  diffraction  and 
Reverse  Monte  Carlo  modeling  allowed  the  amorphous 
structure  of  Ca60MgxZn40_x  alloys  (X=10,  15,  20  and 
25  at.  %)  to  be  studied.  Six  partial  radial  distribution  func¬ 
tions  were  extracted  which  give  a  good  description  of  both 
the  experimental  neutron  and  x-ray  structure  factors  and  total 
correlation  functions  of  these  alloys. 

(2)  Consistent  with  expectations  from  the  efficient  cluster 
packing  model,  the  amorphous  structure  of  Ca-Mg-Zn  MGs 
is  shown  to  be  describable  as  a  mixture  of  Mg-  and  Zn- 
centered  clusters,  with  Ca  dominating  in  the  first  coordina¬ 
tion  shell  of  these  clusters.  A  CN  of  10  [—7  Ca+3(Mg 
+  Zn)]  is  the  most  common  for  Zn-centered  clusters  whereas 
CN=11  and  12  [—7-8  Ca+4(Mg+Zn)]  are  the  most  com¬ 
mon  for  Mg-centered  clusters.  Ca-centered  clusters  typically 
have  CN=  13. 

(3)  Analysis  of  the  Voronoi  polyhedra  in  the  structural 
models  demonstrates  that  there  is  a  link  between  the  clusters 
present  in  crystalline  CaMg2,  and  the  glasses.  As  the  concen¬ 
tration  of  Mg  in  the  glass  is  increased  there  is  a  correspond¬ 
ing  increase  in  the  fraction  of  polyhedra  which  are  both  12 
coordinated  and  have  a  predominance  of  fivefold  bond  ver¬ 
tices.  This  change  drives  the  structure  toward  the  (0,0,12,0) 
polyhedra  found  in  the  crystal  structure  and  explains  the  re¬ 
duction  in  the  GFA.  Similar  changes  are  noticed  in  the  Ca 
environment  but  no  specific  links  can  be  made  between  the 
changes  in  the  Zn  environment  with  composition  and  the 
preferential  formation  of  CaZn. 

(4)  Efficient  atomic  packing  is  shown  around  the  solute 
Mg  and  Zn  atoms,  where  the  Ca,  Mg,  and  Zn  atoms  in  the 
first  coordination  shells  occupy  all  of  the  space  available. 
Packing  around  Ca  atoms  is  less  efficient.  The  degree  of 
packing  efficiency  around  the  solvent  atom  is  suggested  to 
contribute  to  glass-forming  ability. 


(5)  Bond-angle  distributions  suggest  near-equilateral  tri¬ 
angles  and  pentagonal  bipyramids  to  be  the  most  common 
atom  configurations.  While  fivefold  bond  configurations 
(pentagonal  pyramids)  are  still  populous  (—60%)  in  the  first 
coordination  shell  of  the  clusters,  the  characteristic  solute 
(Mg  or  Zn)  centered  clusters  are  not  icosahedral. 

(6)  Similar  to  oxide  glasses  and  a  binary  Pd-Si  BMG, 
Ca-Mg-Zn  glasses  contain  short-range  topological  structures 
similar  to  those  in  their  competing  crystalline  compounds. 
Some  nearest-neighbor  bond  lengths  are  shorter  in  these 
glasses  than  in  the  competing  crystals,  suggesting  that  opti¬ 
mal  short-range  atomic  packing  is  retained  in  metallic 
glasses  since  the  long-range  order  constraint  imposed  by  a 
crystalline  structure  is  not  present. 

(7)  An  ab  initio  MD  simulation,  although  for  a  small 
number  of  atoms,  corroborates  the  structural  findings  above. 
The  extensive  structural  characterization  reported  in  this 
work  thus  offers  a  comprehensive  account  of  the  structure  in 
Ca-Mg-Zn  glasses,  representing  a  category  of  MGs  (alkaline- 
earth-metal  based)  that  has  not  been  well  studied  before. 
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